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We report on a new method to a global fit of the CKM matrix by using the necessary and sufficient 
condition the data have to satisfy in order to find a unitary matrix compatible with them. This 
condition writes as — 1 < cos if < 1 where is the phase that accounts for CP violation. By using 
it we get that the experimental data are to a high degree compatible to unitarity and that ip takes 
values around 90°, in contrast to the previous determinations. Numerical results are provided for 
the CKM matrix entries, the mixing angles between generations and all the angles of the standard 
unitarity triangle. 
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The determination of the Cabibbo-Kobayashi- 
Maskawa matrix that parameterizes the weak charged 
current interactions of quarks is a lively subject in par- 
ticle physics. The four independent parameters of this 
matrix govern all flavor changing transitions of quarks 
in the Standard Model, and their determination is an 
important task for both experimenters and theorists. 
The large interest in the subject is also reflected by the 
workshops organized in the last two years whose main 
subject was the CKM matrix 

Although it seems that at the level of experts there 
exists a consensus concerning the method of extracting 
from experimental data information about the CP violat- 
ing phase 0-01, the usual method of considering only the 
orthogonality of the first and third columns is not reliable 
as we shall show in the paper. The main shortcomings of 
the current approach are the following: a) one uses one 
unitarity triangle instead of six; b) the currently used 
four independent parameters are not re-phasing invari- 
ant; c) one works with an approximation of the standard 
parameterization which may lead to inconsistencies. 

The aim of this paper is to propose an alternative 
method for imposing unitarity by fully exploiting the 
constraints implied by it. We consider that the second- 
generation B-decay experiments of the LHC era, when 
the accuracy will make a tremendous difference, require 
to keep improving all the tools we are working with. Our 
approach provides the necessary and sufficient condition 
the data have to satisfy in order to find a unitary ma- 
trix compatible with them; this condition is given by 
— 1 < cosip < 1, where ip is the phase entering the CKM 
matrix. We construct a theoretical model and we use 
it to test the unitarity property of the data as they are 
provided by the Particle Data Group (PDG) 01 • We find 
that one can reconstruct from PDG data a unitary ma- 
trix and the phase ip is close to 7r/2. 

Our theoretical framework is as follows. We use the 



standard parameterization advocated by PDG that we 
write it in a completely rephaising invariant form as 
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with Cij — COS 6ij and = sin 9ij for the generation 
labels ij ~ 12, 13, 23, and Lp is the phase that encodes 
the breaking of the CP-invariance. It is easily seen that 
by multiplying a unitary matrix at left and at right with 
arbitrary diagonal phase matrices, the phase invariance 
implies that we can freely choose the entries of one col- 
umn and of one row as non- negative numbers. We used 
this property for rewriting the standard parameterization 
in the above form. We remind that phase invariance 
of the J/cATAJ-transformed quark wave functions is a re- 
quirement for physically meaningful quantities. 

On the other hand from experiments one measures a 
matrix whose entries are positive 



(2) 



entries that, in principle, can be determined from the 
weak decays of the relevant quarks, and/or from deep in- 
elastic neutrino scattering [3 . In other words we make a 
clear distinction between the unitary CKM matrix Uckm 
and the positive matrix V provided by the data. The 
main theoretical problem is to sec if from a matrix as 
one can reconstruct a unitary matrix as . If the exper- 
imental data are compatible with unitarity the weakest 
form of this property is expressed as follows 
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We stress that the above relations does not test the 
unitarity, as it is usually stated in many papers; they are 
necessary but not sufhcient conditions. The class of pos- 
itive matrices satisfying Eqs. (O is considerable larger 
than the class of positive matrices coming from unitary 
matrices. The set ^ is known in the mathematical liter- 
ature as doubly stochastic matrices, and the subset com- 
ing from unitary matrices V^^j = \Uij\'^ is known as unis- 
tochastic ones The double stochastic matrices have 
an important property, they are a convex set, i.e. if Vi 
and V2 arc doubly stochastic so is their convex combina- 
tion aVi + {1 — ct)V2, a £ [0, 1] as it is easily checked. 

The first problem to solve is to find a necessary and suf- 
ficient criterion for discrimination between the two sets. 
Assuming the experimental data (|2J) come from a unitary 
matrix as we obtain the following relations between 
the parameters entering Uckm and V matrices 
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the following. The above relations provides the neces- 
sary and sufficient condition the data have to satisfy in 
order the matrix ^ comes from a unitary matrix, and 
this condition is 



1 < cos ip <1 



(6) 
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Since in relations tp enters only in the cosine function 
we can take (ys g [0,7r] without loss of generality. 

The last four relations provide us formulas for cos ip 
and these formulas have to give the same number when 
comparing theory with experiment, by supposing the 
data come from a unitary matrix. Their explicit form 
depends on the independent four parameters we choose 
to parameterize the data, and we will always choose these 
parameters as four experimentally measurable quantities, 
i.e. V^y An other reason is that V'^^ — \U^ckm\^ being 
square of the moduli functions they are rcphaising invari- 
ant; the CP violating phase does not share this prop- 
erty |lflj|-jlll|. By consequence the parameters entering 
on the right hand side of the relation 
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= ■'^13^12^13 + ■S12S23 ~ 2S12S13S23C12C23 coacp 

Vtl = S12S13C23 + C?2S23 + 2.S12S13S23C12C23 COSip 

It is easily seen that the parameterization Q satisfies 
identically the relations J^l, and also that CP-violation 
requires 
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usually used by the physics community o,re not 

rephaising invariant because the left hand side is not. 
Another invariant parameter is the Jarlskog invariant J 
but it it not a measurable quantity. Depending on 
the explicit choice of the four independent parameters 
we get one, two, three or four different expressions for 
cosyj; e.g. if we take Vud, Vus, Vcd, Ks, as independent 
parameters we get 



% ^ 0, ij = 12, 13, 23, 9i2 ^ 7r/2, and ^23 ^ ^/2 



(5) 



We deduce from the factors that multiply cos (p in Eqs.Q 
that an entire region around and 7r/2 is forbidden for 
the above parameters. Using the numerical values from 
PDG data one can get values for cos ip outside the "phys- 
ical" range [—1,1], or even complex as we will show in 
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and from the sixth Eqs.l^J we have 
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In this case, other two independent formulas are given 
by the last two equations in relations Q), and they 
have to give (almost) identical numerical results if the 
data are compatible with unitarity. In general the data 
will give different numerical values for the three func- 



tions expressing cos f. If the independent parameters are 
Vud, Vub, Vcd, Vcb, i.e. we use the information contained in 
the first and the third columns, we obtain four different 
expressions for cos (p, and in this case the mixing angles 
9ij are given by 
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of the four independent parameters. Our proposal is 
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Looking at Eqs. (7) and Q we see that the expressions 
defining the mixing angles are quite different. Thus if 
the data are compatible to the existence of a unitary 
matrix these mixing angles have to be equal, and this is 
a necessary condition for unitarity. As a warning what we 
said before can be summarized as follows: the unitarity 
property is a property of all the CKM matrix elements 
and not the property of a row and/or a column, as it is 
considered by many people working in the field. 

To better understand the above considerations, let us 
consider the toy model defined by = 1/3. It is eas- 
ily seen that no matter how we choose the independent 
variables we get the same numerical results for Sij when 

1/3. In this case we have 



we use the exact values 
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and from any of the last four equations we get cos = 
0, i.e. Lp = 7r/2. Then from Eq.(l) we obtain the unitary 
matrix 
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which is the matrix of the finite Fourier transform in 
d = 3 dimensions. The above matrix has the prop- 
erty that it maximizes the Jarlskog invariant J = 
S12S13S23C12C13C23 sin (p, J„iaa; = l/(6v/3). To see that 
the orthogonality properties of columns and rows are not 
trivial, let us suppose that the second column of the Utoy 
matrix has the form(l, 1, 1)*, where t means transpose. 
It is easily seen that the orthogonality of the first column 
with the second and third ones of this new matrix are ex- 
actly satisfied, but the second and the third columns are 
not orthogonal ! The same happens with all the rows, i.e. 
the orthogonality property is very strong implying sup- 
plementary constraints that until now were never been 
used at their maximum potential. In the real case of 
the experimental data the situation is worse since we use 
approximate values for Vij and then the values of Sij de- 
pend on the choice of the independent parameters that 
define the model and are no more independent of them 
as in the case of the Utoy matrix. 

Now we define a test function that should take into ac- 
count the double stochasticity property expressed by the 
conditions (j2Jl and the fact that in general the numerical 
values of data are such that cos ip depends on the choice 
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Testing our model Hll|l on the toy matrix (|10() we found 
that by adding a second component of the form 



xl-T. E 



(12) 
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the results improve, where Vij is a numerical matrix that 
describes the experimental data, and a is the matrix of 
errors associated to Vij. The test has also shown that 
better results are obtained when the number of cos ip is 
large. In the last sum we use only the data coming from 
the first two rows because the entries of the third row are 
not yet measured. 

The above expressions will be used to test globally the 
unitarity property of the experimental data. 

In the following we test our method on the published 
data, i.e. we want to see if our necessary and sufficient 
criterion, —1 < cosi^ < 1, could constrain enough the 
data. For that we use the PDG data, the fit [3 and 
its recent up-to-date results , and consider that a good 
starting point for a comparison of the unitarity triangle 
approach, almost exclusively used nowadays, with our 
method is to look at the central values from the above 
cited papers. These values arc given in relation l|13|) 
where we used the notation \/V to denote the numerical 
CKM data matrix as it is usually provided. All matrices 
(|13|l satisfy quite well the stochasticity property Q as it 
is seen from the Table 1. 



PDG — 



i/VcifM/G [5] 



CKMfG 



[6] = 



0.97485 0.2225 0.00365 
0.2225 0.974 0.041 
0.009 0.0405 0.99915 

0.97504 0.2221 0.0035 
0.2220 0.97422 0.0408 
0.0079 0.04025 0.99917 

0.97400 0.2265 0.00387 \ 
0.2264 0.97317 0.04113 
0.00826 0.04047 0.999146 / 



(13) 



Here CKMfG denotes the CKM fitter Group. 

As we said before costp depends on the four indepen- 
dent parameters we use to obtain it and for compari- 
son we took six groups of independent parameters in 
the Table 2. It is easily seen that the central values 
from all VpDG and CKMfG are not compatible with 
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TABLE I: ri and Ci, i = 1, 2, 3, are the values on the left 
hand side in Eqs.Q calculated using the data p3^ 

unitarity, although they arc compatible with the double 
stochasticity property to an acceptable accuracy. More, 
the data show that cos ip can take even imaginary val- 
ues and these values can not be properly processed in 
the usual approach of the unitarity triangle and Wolfen- 
stein approximation, which are the starting points of 
all the present analysis of the CKM data. In all the 
cases it is S13 that causes the trouble. For the PDG 
data si3 = + Vtl - ^ud = 4-0 x IQ-^^ i and, re- 

spectively, = ^v^^^+W+v^d + v^ - 10-'*; 

for the CKM fitter Group results [a| it has the form 
si3 = - K?d - ^us = 5-6 X 10-3 i. On the other 
hand such an incompatibility cannot, in principle, be 
detected by using the Wolfenstein parameterization, be- 
cause quantities as S13 = ^ 1 — V^^ — V^^ , that may be- 
come imaginary in some cases, are usually approximated 
by S13 = Vub and never appear in the usual approach. 
The conclusion is that unitarity requires a very fine tun- 
ing between all the entries of the matrix Q and our 
method could put strong constraints on the CKM ma- 
trix. 

For the proper fit we considered six groups of four inde- 
pendent parameters, those appearing in the first column 
of Table 2, that lead to 17 different expressions for cos (p. 
We considered that these six groups could be equiva- 
lent to the six orthogonality relations implied by uni- 
tarity. For minimization of the we used the FindMin- 
imum function provided by Mathematica. Special care 
was taken for properly treating the cases that lead to 
cos if values outside the physical range. 

We started the fit by using the information provided 
by the first group of four parameters in Table 2 since 
we wanted a comparison with the unitarity triangle ap- 
proach that uses the same information. With no indepen- 
dent parameters we found a x' = Xi + X2 = 3.8 x lO"'^ 
and we used the Vij determined parameters to test all 
the seventeen cos(p'*). We found that nine values, for 
i — 1, . . . , 4, 8, 13, . . . , 17, are around 0.05, their mean 
being < cos(p >= 0.0518. However we found also 
a few discrepancies: cos ip'^^^'> — cos ip'^^'' = 0.635 and 
cos lyj*-"^-* = —0.02; even worse four values for i = 9, . . . , 12 
are outside the physical region. We interpret this phe- 
nomenon as showing that the use of only one orthogonal- 
ity constraint leads to non reliable results. 

Taking into account the cosiy9^'' provided by the sec- 
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TABLE II: All the central data ^ coming from 0, [1 and 
are not compatible with unitarity, although the last fit is 
considerably better than the published one. The left column 
contains the six groups of independent parameters used in our 
analysis. 



ond group of independent parameters the fit improves, 
all the cos<^(*^ being inside the physical region, and 
the maximum difference between cosines is of the or- 
der 5. X 10"^ and ip is around 58°, value compatible to 
that provided by the CKM fitter Group |^. Including 
now in all the the first twelve costp'-*-' we found that 
the results considerably improve, and the difference be- 
tween coap^^^ is of the order 2. x 10"^. The big surprise 
was that p takes values close to 90° ! For comparison, 
with our method we obtained by using the up-to-date 
Vij parameters provided by the CKM fitter Group 
that 0.354383 < cos^s^') < 0.558023, values that are 
listed in the last column of Table 2, with a mean value 
< cos(y9 >= 0.460198 and a = 0.0436841 which leads to 
a CP phase p ~ 62.0001°, the interval of variation being 
59.7428° < ip < 65.3853°, results that almost coincide 
with that given in that are 50° < — 7 < 72°. In our 
opinion these results show the limits for the prediction 
power of the unitarity triangle approach. 

If we use only the values provided by the first group in 
Table 2, one gets p ~ 63.7° ±3.°. The difference between 
the two fits comes from their quality. In the following by 
using our method we found that takes values in the 
interval 2. x 10"^ < Xi < 1-7 x 10"^, while the same 
expression provides xl — 2.95 when using the last Vij 
values obtained by CKM fitter Group. 

We included step by step all the constraints implied by 
the six groups and obtained three different matrices; they 
were used to obtain another one by using the convexity 
property of the unistochastic matrices. Other matrices 
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were obtained by relaxing the condition H12|l by taking all 
combinations with five and respectively four Vij that pro- 
vided (^) + (2) = 21 new matrices. The set of these ma- 
trices was considered as 25 independent "experiments" 
on which the statistics was done. For cos (p these lead to 
25 X 17 = 425 values that gave ip = 89.9962° ± 0.0767°. 
The interval for the Xi values for all the 25 matrices was 
shown above and the final results are shown in the Table 
3. In fact the obtained values can be improved. 

One way to do this is to consider the number of all the 
possible four independent parameters Vij which is (^) = 
126, but not all arc independent; we have to exclude all 
groups of four parameters that contain all entries from 
a row or a column, whose number is 36, and also all 
groups that are not independent because of the unitarity, 
e.g. {Vud, Vub, Vcs, Vts), whose number is 9. This leads to 
126 — 36 — 9 = 81 groups and the estimated number of 
cos(p is about 240. The only trouble is that the fitting 
procedure is very time consuming even on a good work 
station. 

Looking at the final results we see that Vij values are 
not far from those provided by other fits. The striking 
feature concerns the values for the angles of the unitarity 
triangle, a, [3 and 7 , angles that are not obtained from 
the fit, but from parameters provided by the fit. One 
sees that with a very good approximation the standard 
unitarity triangle is almost a rectangle one. We have (p > 
7 where their difference is about 4. x 10^^. However this 
leads to the unexpected result, sin2Q! ~ sin 2/3, relation 
that is not satisfied by the experimental data, which in 
particular are not very clean. We obtain that sin 2/3 = 
0.677, value that is not far from the world average of the 
BaBar and Belle experiments, that is sin 2/3 = 
0.736 ± 0.049 Q- In fact both angles a and /3 are at 
the lower extremity of the typical ranges for them , 
70° < a < 130°, 20° < /3 < 30°. We mention that 
a rectangle unitarity triangle was found in |l7l | but with 
a ^ 90°. 

With our values we find Im(U^j^j^,j^^UcKMtd) = 
(1.33 ± 0.03) X 10~^ which is compatible to that used 
in ^3 for determination of the CP-violating ratio e'/e. 

In conclusion we can say that the unitarity triangle 
results are not reliable and that our approach could out- 
perform by far all the other methods used to reconstruct 
a CKM unitary matrix from experimental data. 
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